Higher Dimensional Szekeres' Space-time in 
Brans-Dicke Scalar Tensor Theory 

Asit Banerjee^*, Ujjal Debnath^^ and Subenoy Chakraborty^^ 

^Department of Physics, Jadavpur University, Calcutta-32, India. 
^Department of Mathematics, Jadavpur University, Calcutta-32, India. 
(Dated: February 7, 2008) 

The generalized Szekeres family of solution for quasi-spherical space-time of higher dimensions 
are obtained in the scalar tensor theory of gravitation. Brans-Dicke field equations expressed in 
Dicke's revised units are exhaustively solved for all the subfamilies of the said family. A particular 
group of solutions may also be interpreted as due to the presence of the so-called C-field of Hoyle 
and Narlikar and for a chosen sign of the coupling parameter. The models show either expansion 
from a big bang type of singularity or a collapse with the turning point at a lower bound. There is 
one particular case which starts from the big bang, reaches a maximum and collapses with the in 
course of time to a crunch. 
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The first step to solve the Einstein equations for the metric belonging to Szekeres family [1] 



was taken by Szekeres [1] with dust and A = 0. Here a and (3 are in general functions 
of {t,r,x,y). Szekeres' result was generalized by Szafron et al [2] without any further as- 
sumption except for non-zero pressure. In further generalizations the perfect fluid has been 
replaced by a fluid with heat flow [3], viscosity [4] and also electromagnetic field [5]. Also 
Barrow et al [6] gave solutions for dust model with a cosmological constant and recently 
Chakraborty et al [7] gave solutions for perfect fluid model with a cosmological constant in 
(n-|-2)-D space-time. We are not aware of any generalization of the above class in more than 
4 dimensions in Brans-Dicke scalar tensor theory of gravitation. In this paper, we work out 
solutions for dust in the presence of the cosmological constant (A ^ 0) and the Brans-Dicke 
scalar field. We consider the scalar tensor theory in the Dicke's [8] revised version after unit 
transformation. In the revised version, G does not vary while the masses of the elementary 
particles are varying. In this version, the trajectories of particles are not geodesies and the 
energy momentum tensor holds for combined matter and the scalar field. Suitable transfor- 
mation of units used by Dicke are (jf^u = ^g^iv, ^ — A~^/^to, ds = X^^^ds and the scalar 
field A = (/>/</'0) where bars indicate the variables in the revised units, A is the scalar field in 
the new unit and (f>o is a constant. The field equations in the revised units are 



where S^j^i, = ^^g^g^ {Cf^Cu — ^g^^^CiC^) and the scalar C = In A = In (0/(/)o). 

In what follows we consider a more general (n + 2)-dimensional Szekeres space-time, which 
can be expressed by the following metric form [7] 



I. INTRODUCTION 



ds^ - dt^ - e'^dr^ - e^^idx^ + dy^) 



(1) 



(2) 



n 




(3) 
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where a and (3 are now functions of all the (n + 2) space-time variables i.e., 

a = a{t,r,xi, ...,a;„), (3 = (3{t,r,xi, ...,a;„). 

It is to be noted here that in the presence of the non-zero scalar field it is necessary to 
make at least one assumption in order to solve the field equations. We can either assume 
the form of scalar field or of the metric co-efficients. Ifere we have assumed =0. Then 
there arc two possibilities. In one of these cases the scalar field is in general a function of r 
and t, whereas in the second case the scalar field may either be a function of r or a function 
of t alone. In fact the differential equation obtained in terms of the metric co-efficients for 
C = C{t,r) is apparently not solvable. So in the next section we present exact solutions in 
the following different cases: 

(i) =0,C = C(t,r), 

(ii) /3V0,4, = 0,C = C(t), 

(iii) p' = 0,p,,=0,C = C{t). 



The case C = C{r) is not consistent with the field equations. In all the above cases the 
asymptotic behaviour of the dust density p and the shear scalar a are obtained in the limits. 
The cosmological models described in all the cases mentioned above are all inhomogeneous 
models. 



The special cases of our solutions in the absence of the scalar field lead to those of 

Szekeres as shown in the subsequent sections. At this point one must note that Szck- 
eres solutions are however given in co-ordinates different from those used in the present text. 

In the last section an alternative interpretation is given for the above mentioned solutions 
in the special case of C as a linear function of time. In this particular case the scalar field 
C may be interpreted as the creation field first proposed by Hoyle and Narlikar [9] in order 
to obtain a steady state cosmology or later a quasi steady state cosmology. 



II. FIELD EQUATIONS IN THE REVISED VERSION OF THE 
BRANS-DICKE THEORY AND THEIR EXACT SOLUTIONS 

The Einstein's field equations in the presence of the cosmological constant A follow from 
(2) (we choose SwG = 1, see [7]): 

1 " f 1 

nd/3 + -n(n - - e"'^ | + 2 ^" " ^^^"^ ~ ^^-^'^ + ~ 2)a,,/3,, + a,,,, 

+(n - l)/J,,,J + e-2« La' (3' - Jn(n + l)/3'2 - n(3"] = A + p- i/ (c^ + e'^^C'^) (4) 



1 1 " r 1 1 

-n(n + +nP- -n{n - l)e-'"P'' - e'^^' E o ~ ^^^"^ ~ ^^^'^ + ~ 

= A+i/(c7^ + e-2«C'2) (5) 



op- +a^{n- -n(n - l)/?^ (n - I)/3 e^^a j _ Y)^^ ^ _ _ 
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(n - 1)/?"} - e-^^ {< + - 2)(« - ml, + ax,x, + (n - 2)/3,,,, + (n - 3)a,,/J,, | 
|(n - l)a,,/3,, + i(n - l)(n - 2)/32^| = K+\f {c^ + e-'"C'') (6) 

Q!x,- (-axi + Pxi) + Pxj (ttx, + {n- 2)/3^,) - a^.^. - (n - 2)P^,^. =0, (i 7^ j) (7) 



-2/3 



a/3'-/3/3'-/3' = -l/C'C7' (8) 



- aa^, + Pa^, - a^, - (n - 1)^^, = (9) 



=0 (10) 



where dot, dash and subscript stands for partial differentiation with respect to t, r 
and the corresponding variables respectively (e.g., (ix^ = ^) with z,j = l,2,...,n and 
/ = -i(2a; + 3). 



From equations (8) and (10) after differentiating with respect to Xi and t respectively, we 
have the integrability condition 

iix,0^ = ~fCC'p'x,, i = l,2,...,n (11) 



This equation cannot be solved without any specific assumption either on the metric 
co-efEcients or on the scalar field C. Such an assumption is however necessary to obtain 
exact solutions for all the variables. One of such assumption is /3^. = 0, which leave us two 
possibilities: (i) 7^ so that (ixi = 0, (ii) /?' = 0. The second choice in view of eq.(8) 
leads us to the condition CC = 0, which implies that C can not be a function of both the 
co-ordinates r and t. For general C ~ C{t,r), first of all we consider /?' 7^ 0,$xi = 0. So 
we obtain from the field equations (8) and (10), the metric co-efHcients are in the following 
form: 



e'' = i?(t,r) e'^'^'^i'-'^") (12) 

and 

R' + Rv' 



e 



S{t,r) 



(13) 



where R and S arc function of t, r only. Now from the field equations (5) and (6) using 
equations (12) and (13), we have the differential equations for R and S: 

2RR +(n- l)(i?2 _ 62^2) _ IpA + f(C^ + e-^"C'^)]R'' = (n - l)K(r) (14) 

n 



4 



and 



2n 



(15) 



where JsT is a function of r only. Since (3'^. =0, so a is a function of t, r only and the 
form of e'^ suggests 



(16) 



where A{r) is arbitrary function of r alone satisfying {n — 4)j4^(r) = K{r). It is to be 
noted that when n ^ A then K{r) ^ and n = 4 implies K{r) = 0, which shows that 
K{r) must vanish when we consider a six dimensional space-time. In other cases, however 
K{r) ^ 0. The field equation (7) is automatically satisfies by the above solutions. Also from 
the field equation (4), the expression for the energy density p is 



p{t, r, XI, x„) = ^(a + nl3 + a' + np^) + + -^{C^ + e'^^C'^) (17) 

n — 1 „ T „ T 



n— 1 n — 1 



and the expression for the shear scalar is 



{RR' - RR') 



l\2 



2{n + l)R'^{R' + Rv'Y 



(18) 



Now eliminating the terms containing the derivative of the scalar field C between (14) 
and (15) we obtain the differential equation in R and S as 



d_ 

dr 



[2RR+ (n - - - AT)} 



This differential equation in R and S is quite complicated and is apparently not solvable 
in closed form. For general C = C{t,r), it is not possible to find out any exact solutions. 
So wc may consider C is a fimction of t or r alone. For C = C{t), wc arrive the following 
two cases from equation (11) in order to obtain explicit form of the metric co-efficients: (i) 
/?' ^ 0, = 0, (ii) /?' = 0. 



Case I: /?' 7^ 0, = (i = 1, 2, n), C = C{t) : 

With the above condition being used we obtain from the field equations (8) and (10), the 
metric coefficients in the following form [7]: 



e'^ = R{t, r) e 



i/(r,a;i,."iiCn) 



(20) 



and 



e" = R' + Ri^' 



(21) 



Now using the relation (20) and (21) in equation (5), we obtain a differential equation for 
R as follows 



2RR +(n- 1)R^ --(2A + fC^)R^ = {n- l)K (22) 
n 



and the solution for v can be expressed as 



e-'^ = ^^x2 + ^SiXi + L» . (23) 



j=l i=l 

There is, however, the following restriction 



^ B} - 4:AD = K-1 (24) 



where K, A, Bj's and D are all functions of the co-ordinate 'r' alone. 
The first integral of the equation (22) may be given by 

P OA 1 r 

E? = K+ -— + + - / C'^RJ'Rdt (25) 

n{n + l) n'' J ^ ' 

where F is an arbitrary functions of 'r' and / is a constant which is already expressed 
earlier in terms of Brans-Dicke parameter w. 

The following situations are considered separately: 

(i) C = Cot, K ^ and n = 3, that is this satisfied A = e*-^"* where Cq is a constant. 
The integration of the equation (25) in the next step yields 

2zR^ + K = {AzF - K"^)^^^ Sinh{2z^^^{t-to)}, z>0 (26) 

and 

2\z\R^-K={K^+4\z\F)^/^ Sin{2\z\^/^{t-to)}, z<0 (27) 
In the above solution for R is another arbitrary function of r and is a constant given 

•' n(n+l) 

(ii) C = Cot, K = 0, n is arbitrary: 
The solution is given by 

1 f F\^^'^ 1 
R^=i-\ Sinh{-{n+l)z^''^{t-to)}, z>Q (28) 



and 



1 /2 

^"^"(l^) Sin{\{n + l)\z\^l\t-to)}, z<0 (29) 

Here of course there is another possibility. If F < 0, we can obtain from (25) on integration 
a different solution expressed as 



n. + l 



F\V'^ ^ . 1 



= I^^J Cosh{^{n + l)z^'^{t-to)}, z>Q (30) 
The solutions (28) and (30) have different properties and will be discussed subsequently. 
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(iii) C = Cot", K = 0, 'a' being a constant parameter. We put R = G^/("+^) in the 
equation (22) so that we obtain 

AnS-^G - (n + 1)(2A + fC^) = 
which on assume A = may also be written as 



G - yr"--^G = (31) 

where y = ^"4^^^ fC^a?- In equation (31) after integration finally yields the following 
solution for the variable R, 



G = R"^ = Ui(r) /i[yi/'r/a] +^2(r) I.^iy^^H^/a] 



(32) 



and 



G = i?^ = t'/^ U (r) J^[\y\'/Hya] + A2{r) J^[\y\'/Hya]\ , y<0 (33) 



It may be noted that when a = 1, the solutions (32) reduces to either (28) or (30) but 
solution (33) reduces to (29) only. As usual Jm{x) and Im{x) respectively stand for Bessel 
function and modified Bessel function of first kind of order m. 



The expressions for the density p and the shear scalar a are obtained from (17) and (18) 
and are given by 



p{t,r,xi,.. 



n R' + Rv' nR 
(n - 1) R' + Rv'^^ 



+ 



(2A + nfC^ 
(n-1) 



(34) 



and 



2 ^ n {RR! - RR'f 
" 2{n + l)R^{R' + Ru')^ ^ ' 

The explicit forms of (34) and (35) arc apparently very complicated. However one can 
discuss the behaviour of different models in the limits t ^ to and t ^ 00. In fact to{r) 
depends the initial moment of evolution in each case. The evolution may or may not begin 
with a big bang singularity. Again since 7^ the instant of singularity if it exists will be 
position dependent. These are apparent from the following facts. 

In case (i) with K ^ the dust density and the shear scalar both are finite at t ^ to but 
in the other extreme limit that is as t ^ 00 the dust density p remains finite but the shear 
scalar a vanishes. 



In case (ii), p ^ 00 and cr ^ 00 at the initial instant t ^ to but in course of time as 
t ^ 00 the dust density remains finite, whereas ci ^ 0. This is true for the solution (28). 
On the other hand if we concentrate our attention on the solution (30). We find that the 
density p remains finite throughout the evolution, whereas the shear starts from a finite 
magnitude and gradually disappear in course of evolution. 

In case (iii), we note that there is singularity {R = 0) for the solution in equations (32) 
and (33) provided we choose arbitrary constants A2{r) = 0. Also asymptotically for large 
t, R oscillates infinitely for the solution (33) while R becomes infinite for the solution 
(32). In this case, for < a < 1, p and a both explode initially whereas they attain finite 
magnitudes in course of evolution as t — > 00. 



7 



Case II: /?' = 0, C = C{t): 

In this case, for only the above choices, the exact solutions can not be found from 
the field equations (that shown by Szekeres [1]). So we need further assumption like 
$xi = (i = 1, 2, ...,n), in order to obtain the exact solutions. 

Now from the field equations we obtain the metric coefficients in the form 



e'^ = R{t) e' 



1^(0:1, X2,.">^n) 



(36) 



and 



e" = i?(t) T]{r,Xl,X2,...,Xn)+^iit,r) (37) 
Then as before from the field equation (5), we have similar differential equations in R as 



1 



2RR +{n- 1)R' - -(2A + /C^)i?^ = (n - 1)K 



and the solution for u as 



along with the restriction 



(38) 



(39) 



^ Bf - AAD = K 

i=l 

Here A, Bi's and D are all arbitrary constants. 

Now from the field equation (7) we have the solution for rj as 



(40) 



i=l i=l 

and the resulting diff'erential equation in fi is 



(41) 



with 



R'jl + (n — l)i?/i + ji 



g(r) = (n - 1) 



1 



R-- (2A+/C^)i? 
n 



2{uD + wA) -^ViBi 

i=l 



9{r) 



(42) 



(43) 



where u, ViS and w as arbitrary functions of the co-ordinate 'r' alone. 

For simplicity, let us choose C = Cot, n = 3. In this case the solutions of R (as before) 
and n are 



22i?2 + K= {4zF - K'^f'^ Sinh{2z^''^{t - to)}, z > Q, 



(44) 
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2\z\R'^-K={K^-AzFy/^ Sin{2\z\^^^{t-to)}, z<0 (45) 

and 

2ziJLR + g{r) = ^/Az'^h{r) - g'^{r) Sink |\/2i [t - ii(r))| , > 0, (46) 

2\z\^iR-g{r) = ^/ g'^{r) - Az'^h{r) Sin\y2\z\ (i-ti(r))| , ^ < 0, (47) 

where Iq, F are arbitrary constants, h{r), ti{r) arc arbitrary functions of 'r' only and 
z = ^(2A + fC^). The solutions (44) and (45) is identical with (28) and (29) except for 
the fact that now F and to are no longer functions of r. 

In this case the expressions for the density and the shear scalar are given by 



p{t,r,xi, ...,Xn) 



(n-1) 



il + Rt] 
IJ. + Rr] 



nR 



(2A + n/C^) 
in -I) 



(48) 



{R12 - Riif 



2{n + l)R^{n + Rr])^ ^^^^ 

There arc two distinct cases for cither K or K = 0. In the former case with K < 
both the dust density and the shear scalar remain finite when t ^ to, but they are infinitely 
large at this instant in the second case = 0. On the other hand as t ^ 00 the density 
is finite even though the shear vanishes in course of evolution. When K > the situation 
differs. Here the density and shear explode at same initial instant other than t — to. 



Next if we consider C = C{r) then from equation (11), we may get also two possibilities: 
(i) P' 0, — 0, (ii) /?' = 0. In these cases the form of metric co-efficients arc similar to 
case I and Case II, but the differential equations (22) and (38) are slightly different i.e., 

2RR +{n- l)R? - ^(2A + fe-^"C'^)R^ = {n- 1)K 

This equation can not be solved because q is a function of all space-time co-ordinates. 
So we are not going to further discussion on the choice C = C(r). 



III. ALTERNATIVE INTERPRETATION OF THE SOLUTIONS IN 
TERMS OF C-FIELD COSMOLOGY 

It is appropriate at this stage to interpret some of the previously described solutions as 
due to the presence of the creation fi(^ld first introduced by Hoylc and Narlikar [9] . The field 
equations in this case are exactly identical with (2) except for the replacement of S^i, by 
Tji'i), where 

T^^'J = -/ (c^C. - i g^^CiC'^ (50) 

where = C,^ and / (> 0) is a coupling constant. The additional feature in C-field 
cosmology is that one must confirm that the C-field satisfies the source equation 

(51) 

It is not difficult to check that the Bianchi identity and the source equation (51) lead to 
the relation C — 1, which determines the expression of the C-field scalar C — t + tp{r). So 
for obvious reasons all the previous solutions with the scalar field C = In A expressed as a 
linear function of time are also solutions of the C-field cosmology. 
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IV. DISCUSSIONS 

One must note that in all the models presents above some are singularity free and some 
have big bang type singularities at the starting point. Particularly for the function K = 
those which show the turning point {R = 0) at some stage represent only the lower bound 
(since i? > 0) as are evident from from the equations (22) and (38). The singularity 
occurs in each case at an instant which is not fixed for different shells of different r. In 
fact, the singularity i? = is position dependent in all the cases described above and the 
reason is that the system is inhomogeneous which is the limit of spherical dust reduces to 
Tolman-Bondi space-time. There is one particular case given by the solution (29) where 
there is an initial singularity at i — tg. Subsequently R{t,r) increases and reaches a 
maximum followed by a collapse to a crunch. 

Another interesting point to observe is that for a few solutions such as (26), (28), (44) 
etc, we have R initially less than zero but becomes positive in course of time indicating 
that the expansion starts from decelerating phase to an accelerating phase at late stage. 
It is to be noted that in the above set of solutions is a kind of singularity given by 
e" = or equivalently (e'^) = 0, which are analogous to the shell crossing singularity 
in Tolman-Bondi [10] models. In this case also the density p diverges. Finally it can be 
mentioned that one can obtain Szekeres 4 dimensional solutions from our solutions if we 
put n = 2,A = and C = 0. In fact the differential equations in R (see equations (22) 
and (38)) after the above simplification becomes identical to the corresponding differential 
equations in Szekeres solutions [1]. Therefore our solutions are generalization of Szekeres 
solutions. 
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